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Abstract 

By repeated application of the deletion-contraction identity we compute the partition function of the Potts 
model on lattices with recursive symmetry with arbitrary values of q and temperature parameter v — e K — 1. 
We first apply the technique to analyze strips of width L y = 2, for three different lattices: square, kagome and 
'shortest-path' (to be defined in the text). It is shown that this approach is not restricted to the explicitly 
solved lattices, but can be applied to a very broad family of graphs of finite width constructed with an arbitrary 
number of recursive steps. Additionally, we show that the same procedure can be applied to wider strips by 
solving the kagome strip of width L y — 5 with the aid of a computer program. As a demonstration of the 
versatility of the method we construct the exact solution of a recursive lattice similar to strips but whose width 
changes periodically from L y = 3 to L y — m, where m is an arbitrary integer number. Lastly, we analyze a 
non-periodic lattice for which the transfer matrix becomes step-dependent. 

Keyword: Potts model, exact results, kagome lattice, diced lattice, dichromatic polynomial. 
PACS: 12.40.Nn,11.55.Jy, 05.20.-y, 05.70.Fh. 

1 Introduction 

The Potts model [T] is one of the most interesting models in statistical mechanics. It appears as the most natural 
generalization of the Ising model by allowing the spin variable to take more than two distinct values. Besides its 
intrinsic importance in the theory of critical phenomena, it manifests a number of intriguing relations with many 
areas of both physics and mathematics (for nice detailed reviews see [21 [31 HI El E] ; for the relation with the problem 
of color confinement see [7] ; for its relation with the Kohanov Homology see [5] ) . 

We will study the q— state Potts model at zero magnetic field, whose partition function on an arbitrary graph 
G is, 

Z(G,g,/?) = 5>-^ 

where /? = 1/ksT and 

m 

Spin variables <Ji can take values 1, q and are located at the vertices of G. There is also an interaction energy J 
for pairs of spins located at the ends of every edge (ij) G G. It will be important for our development to introduce 
an equivalent expression for the partition function, 

Z(G,q,v)= £n( 1 + t,< Wi)' ( L1 ) 

where the temperature variable v = e K — 1 has been introduced with K = f3J. For the ferromagnetic case (J > 0) 
one has < v < oo corresponding to the temperature interval oo > T > 0; for the antiferromagnetic Potts (J < 0) 
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the interval — 1 < v < corresponds to < T < oo. It can be shown directly from (jl.ll) that the partition function 
admits also the following polynomial (Fortuin-Kasteleyn) representation [9] [10] , 

Z(G,q,v) = J2 q^\<^ (1.2) 

G'CG 

where G' is a graph with the same vertex set V of G and an edge set E' C £7, with £7 being the edge set of G; fc(G') 
and e(G') are respectively the number of connected components and the number of vertices of G'. It is important 
to notice that using (|1.2[) one can easily extend q and v from their physical values (q € Z + , v G [—1, +oo)) to the 
whole complex plane. 

Unlike the case of the Ising model (q = 2), solved by Onsager [11], the exact solution for the Potts model in any 
infinite two dimensional lattice is still not available. Nevertheless, thanks to universality, the critical behavior of 
the ferromagnetic Potts model is fairly well understood. On the other hand, the antifcrromagnetic Potts depends 
on the structural properties of the lattice and its critical properties can vary strongly from case to case. Thus, the 
search for tools to obtain analytical information on the free energy of the Potts model on generic lattices is a very 
important task. 

An approach that has proven fruitful in recent years is to solve the simpler problem of calculating the partition 
function for periodic strips of finite width (L y ) and arbitrary length (L x ). Previous work along these lines has 
relied basically on a transfer matrix method [HI [T31 HH HH US] to obtain exact results at arbitrary temperatures 
for strips of widths going up to L y = 7 (honeycomb lattice with free boundary conditions) [17] • Such efforts have 
produced detailed studies for the square, triangular, and honeycomb lattices [P8 ] IT4 ] [15 ] [T6 ] [19 | [20 ] IT? ] I2T]. 

In many cases it is useful to exploit the recursive structure which is present in lattices of physical interest. Even 
when the analytic solution is not available, using a simple ansatz which respects the recursive symmetry one can 
get an excellent agreement with the numerical data [35] [23] [55] . Motivated by these facts we develop a method 
that permits the calculation of the exact partition function for strips of virtually any kind of lattice. 

The idea to use the dichromatic polynomial in the cases of recursive lattices has been proposed in (26] [27] [28] , 
where the authors derive a scalar homogeneous recurrence equation of order higher than one, whose solution is the 
dichromatic polynomial of the corresponding graph. In the present paper it will be shown that it is possible to 
extend this idea by formulating vectorial homogeneous order one recurrence equations whose solutions represent 
the dichromatic polynomials of the graph corresponding to different boundary conditions. As it will be shown in 
the following sections, this vectorial formulation also allows to deal with recursive lattices which are less regular 
than strips and whose shapes change with the recursive steps. 

The paper is organized as follows: In Sec. [2] we explain the general method and illustrate its use by solving 
the L y = 2 and L y = 3 square ladders for which the results are already known [T3l [16]. Next, in Sec. Owe 
present new exact results for arbitrarily long strips of width L y = 2 for the kagome, diced and 'shortest path0 
lattices. In order to demonstrate the versatility of the procedure, in Sec. [5] we obtain the exact solution of a L y = 3 
strip corresponding to a mixture of triangular and square lattices, which additionally has periodic attachments of 
segments of a different strip. Lastly, in Sec. [6] we analyze an interesting case of a non-periodic lattice for which the 
transfer matrix becomes step-dependent. Conclusions and prospects for future work are discussed in Sec. [7] 



2 Recursive equation method 

The Potts model partition function on a graph G (Z G ) satisfies the well known deletion-contraction theorem, which 
says that it can be written in terms of Z Gd and Z Ga as, 

Z G = Z Gi + vZ Gc , (2.1) 

where Gd corresponds to the resulting graph when deleting an arbitrarily chosen edge from G; and G c is ob- 
tained from G by contracting the same edge (i.e. deleting the edge and identifying the vertices at the end of it). 
Furthermore, it is easy to show that the Potts model partition function also fulfills the following property, 

Z { . G ) = qZ Gl (2.2) 

where Z^ m G ) denotes the partition function of a disjoint union of an isolated vertex and a graph G. The resulting 
partition function is a polynomial in q and v which in the the mathematical literature is known as the dichromatic 
polynomial of G. 

lr This lattice will be introduced in Sec. 13.31 It corresponds to a square lattice plus a graph inside each square that has a topology 
identical to that of the shortest path joining all four corners. 
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Figure 1: Basic building blocks for a periodic lattice strip, (a) H corresponds to the structure that will be repeated 
n times. Empty circles to the right denote vertices that will be shared between adjacent blocks, (b) Graph I at 
the right end of the strip containing the information of the initial condition Z(0) = Zi. (c) The strip is built up by 
repeatedly joining ii-blocks after /. 



Identities (|2.1[) and (|2.2| will form the basis of the method developed in what follows. 

Using the notation of [29[[30], consider now strip graphs of the form (G s ) n = (11™= 1 where H is repeatedly 
attached after the initial subgraph /, which will be put to the right by convention. The junction between contiguous 
subgraphs is done by sharing a subset of L y vertices (V") as illustrated in Fig. [TJ Each block / or H can contain an 
arbitrary number of vertices besides the ones in V. Suppose now that starting from a strip of length n as the one 
shown in Fig. [TJ:, we apply repeatedly the deletion-contraction theorem to all the edges in the last H block. Once 
this process is finished, it is easy to sec that the partition function of the initial graph (Zo(n)) will be written as a 
linear combination of the partition functions of objects of length n — 1, 

N-l 

z o(n) = aoj(q,v)Zj(n-l). 

3=0 

The coefficients aj(q, v) are polynomials in q and v and appear naturally as the edges are contracted (factor v) and 
vertices become completely disconnected (factor q). A number N of graphs of length n — 1 result once all of the 
edges in the n-th £f-block are deleted or contracted. Evidently, these configurations arise from the different ways 
in which the L y vertices in V between the n-th and (n — l)-th i7-blocks become identified with each other, by the 
action of the deletion-contraction theorem over the n-th block. One can proceed analogously starting with each 
one of the possible Zj(n); the emerging linear system of equations can be written as 

Z(n)=A(q,v)Z(n-l), (2.3) 

where A(q,v) is an N by N matrix with elements a,ij(q,v), and the vector Z(n) arrange the partition functions 
for all the various possible connectivities at the top layer. It is worth note that A(g, v) is independent of the step 
n: this is a great simplification related to the particular recursive symmetry of the system. As it will be explained 
in detail below, the dimension TV of matrix can be calculated considering only some basic properties of H (i.e. 
symmetries or if it is planar or not). By recursive use of (|2.3p . it is straightforward to show that the general solution 
is 

Z(n)=A n (q,v)Z(0), (2.4) 

where Z(0) is the partition function vector corresponding to the step— block /, in the context of the recursive 
equations (|2.3|) Z(0) play the roll of initial conditions. Wc can now use standard linear algebra methods to rewrite 
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H = 



(a) 



I = 



(b) 



n n-1 2 1 

(c) 

Figure 2: For the square ladder it is quite simple to find its basic building blocks, (a) Every H adds two new 
vertices (full circles) and three edges. A block is joined to the previous one by sharing the L y = 2 vertices to the 
left (empty circles), (b) In this case the initial graph I consists of two joined vertices. 

the solution (|2.4p in the form, 

s=0 

where A s (a s ) are the eigenvalues (eigenvectors) of A, and c s are some constants sensible to the vector of 0— step 
lattices Z(0) = Zj. 

Let us now summarize the method by reviewing its basic steps. 

Step 1: Identify the periodic structure of the strip and choose basic units (H) which are joined with each other 
by a set of L y vertices. 

Step 2: Using the deletion-contraction theorem, write down the partition function of a strip of length n as a 
linear combination of graphs with length n — 1. These new graphs will correspond to different ways in which the 
vertices at the top layer become identified with each other as a result of the deletion-contraction procedure. 

Step 3: Repeat Step 2 for all of the possible connectivities at the top layer. The result will be a matrix recurrence 
equation whose solution in the form (|2.4|) contains the desired partition function. 

This approach has been proposed and applied for zero temperature calculations in [3TJ [301 HI]- Also, it is 
important to notice that to find A(g, v) is equivalent to construct the transfer matrix in the basis of connectivities 
at the top layer as explained in [14] . Nevertheless, the conventional transfer matrix approach requires a layer by 
layer construction of horizontal and vertical matrices that strongly restricts the manageable topologies. A clear 
advantage of the framework presented here, is that it becomes easier to handle complicated lattice topologies, as it 
will be illustrated below. Lastly, it is worth noticing that the algorithm developed above can be implemented in a 
computer software in order to handle larger blockso 



2.1 Application to L y = 2 square ladder 

Here we will consider the case of 2 x n square ladder latticesjf] As illustrated in Fig. [21 it is easy to find for this 
lattice its corresponding initial graph / and its building block For this case we have L y = 2 and by direct 
calculation it is easy to see that N = 2: as explained in Fig. [3J we only have here the original graph corresponding 
to Zq and one graph with the two vertices at the end identified with each other. The vector of the initial graph / 
is then, 

z » =(:!;::!)• <«> 

In order to construct matrix A sq (g, v) for the square (sq) strip, we need to delete and contract all three edges in 
H as it is explicitly shown and explained in Fig. 01 A similar procedure is followed starting from Z\(n) to obtain 

2 In fact, the authors have already developed such a code for L y = 2 systems, in which the user only needs to input the desired graph 
for the building bock H and receives the transfer matrix as an output. Extensions to solve wider strips are obviously possible. 
3 We solve this system for illustration purposes only since the solution was already obtained in |13| . 

4 It is important to observe that there is a certain arbitrariness in the choice of H. For example, one could choose instead to join the 
graphs by the two vertices along a diagonal of the square. 
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Z„(n) 



Z,(n) 



(a) 



n 

(b) 



n-1 



Figure 3: (a) Za(n) is the sought partition function which in our convention is the first component of Z(n). None 
of the end vertices of the graph are identified with each other, (b) Once the deletion-contraction theorem is applied 
to the previous block, a new kind of graph will appear in which the two leftmost vertices arc identified with each 
bther. In our notation identified vertices are joined by a crossed line. 
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Figure 4: The deletion-contraction theorem (|2.1[) is applied successively to each edge in the last block. As a result, 
the partition function of the original graph of length n is written as a linear combination of the partition functions 
of two distinctive objects of length n — 1. 



finally, 



(q + v)(q + 2v) + v 2 



(l + v)(q + 2v) (l+v)v 2 
It is now trivial to diagonalize this 2x2 matrix to find that its eigenvalues are 

A s ± q (<7,v) = i(q 2 + 3qv + 4v 2 + v 3 ± vV + 6q 3 v + 13q 2 v 2 + 16qv 3 - 2q 2 v 3 + 12w 4 - 2qv i + Av 5 + v e ). (2.6) 

As it is well known, the dominant eigenvalue A + will determine the free energy in the thermodynamic limit. It can 
be noted that the discriminant is always positive in the physical regime. 



2.2 Solution of the L y = 3 square lattice 

In order to illustrate how our method can be also used for lattices with L y > 2 we solve here explicitly the 3 x n 
square lattice. The exact solution for such a system was first derived in [16] . 

The idea, as before, is to use rules (|2.1|) and ([2.2)1 to derive a relation between the partition functions of the 
lattices of length n and n — 1: 

^(n)=A^(g,«)^ ) (n-l), (2.7) 



where, 



A (3)(<7> W ) 



I an 012 aw ^ 

0,21 a-12 «23 a 24 
«31 a 32 ^33 a 34 

\ a 4 i a 42 043 a 4 4 J 

The four components of ^(3) are the partition functions on a 3 x n square ladder lattice in which on the rightmost 
vertices the four basic identification have been respectively performed (see Fig. [5]). As in the previous subsection, 
the functions are easily computed using the rules ()2.ip and ([2.2p . The exact values of the functions are 
indicated in Table [TJ 
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Z7(n) 



Z7(n) 



Z'I(n) 



Z" s '(n) 



Figure 5: The four components of the partition function vector correspond to the different ways the vertices at the 
end of the lattice can be identified with each other as a result of the deletion contraction procedure. 

a n = (q + 2v)(q 2 + 3qv + 4v 2 ), a 12 = 2v 3 (q + 2v), a 13 = v 5 , a 14 =v 4 , 

a 21 = (l + v){q 2 +4qv + 5v 2 ), a 22 = v 2 (l + v){q + 3v), a 2S = v 4 (1 + v) , a 24 = v [i (l + v), 

a 31 = [q + 3v)(l + v) 2 , a 32 = 2v 2 {1 + v) 2 , a 33 = v 3 (l + v) 2 , a 34 = v 2 (1 + v) 2 , 

an = q 2 + 5qv + 8v 2 + qv 2 + 3t> 3 , a i2 = 2v 3 (2 + v) , a 43 = v 4 (2 + v) , a 4 4 = v 2 (q + 3v + v 2 ). 

Table 1: Values of the functions ay. 



One can easily obtain the initial condition by building the vector ^^(0) of the partition functions of the 
elementary 3x1 ladders shown in Fig. [51 



%0) = 



( q(q + v) 2 \ 

q{q + v)(l+v) 
q{l + v) 2 
\ q(q + 2v + v 2 ) J 



(2- 



The resulting partition function for any n is then: 

^(31 ( n ) 



A sq 

A (3) 



2(2) (°) 



Thus, the present method provides one with results which are equivalent to the ones obtained with more usual 
methods [15] . 

3 Some new exact results for L y = 2 strips 

In order to illustrate the power of the method described above, we will use it now to obtain the exact partition 
function for strips of three different lattices: kagome, diced and 'shortest-path' (described below). To analyze the 
results we plot for each case the internal energy per site and specific heat as a function of the temperature variable 
v, for different values of q. 

Internal energy and specific heat are quantities that can be obtained straightforwardly from the free energy by 
taking the appropriate derivatives as follows: 



TP df T( ^U d f 



and 



K 2 {v + 1) 



C _ 1 dE 
Kb Kb dT 

In the limit of infinite length the free energy per site becomes 



df | | d 2 p 
dv dv 2 



1 



f = In A 



(3.1) 



(3.2) 



(3.3) 



where N uc is the number of sites per unit cell and Ad is the dominant eigenvalue of the transfer matrix. For 
simplicity, instead of the internal energy, we will plot the dimensionless quantity —E/J, which has the virtue of 
having the same sign for both the ferromagnetic (0 < v < oo) and antiferromagnetic (—1 < v < 0) physical regions. 



3.1 Kagome lattice strip 

The kagome lattice is the medial graph of both the triangular and hexagonal lattices. The Potts model on this 
network is a prototypical highly frustrated system and has been widely studied for many years. Nevertheless, in 
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Figure 6: These partitions functions of these four graphs constitute the initial condition vector Z^(0) 





Figure 7: A strip of the kagome lattice with L y 
The transfer matrix for this unit cell (eq. (|3. 410 
theorem. 



= 2. Each unit cell is composed of five vertices and eight edges, 
is obtained by successive application of the deletion-contraction 



contrast with simpler lattices like the square, triangular or honeycomb [32J [33j [34] , many fundamental questions 
such as the exact determination of the critical frontier remain unresolved in the kagome network [351 133 [33 I3H] • 
In this context, exact results for kagome strips of finite width may provide some insight into the properties of the 
model in two dimensions. 

We proceed to use our recursive equation method to calculate the exact partition function for the L y = 2 kagome 
strip shown in Fig. [7] The resulting transfer matrix is, 



where, 







0-12 \ 


(3.4) 




\ "21 


4a )' 


k 
11 


= q 5 + 8q 4 v + 28q 3 v 2 + 54gV + 2q 3 v 3 + 59qv 4 


+ 10q 2 v 4 + 30v 5 + 20qv 5 + 16v 6 + qv 6 + 2v 7 , 


(3.5) 


k 
12 


= q 2 v 4 + 6qv 5 + 9v 6 + 2qv 6 + 6v 7 + v a , 




(3.6) 


k 
21 


= q 4 + 8q 3 v + 27q 2 v 2 + q 3 v 2 + 46qv 3 + 10q 2 v 3 -+ 


36t; 4 + 30qv 4 + 2q 2 v 4 


(3.7) 




+34,v 5 + lOqv 5 + Uv 6 + qv 6 + 2v 7 , 




(3.8) 


k 
22 


= 2qv 4 + I2v 5 + 13v 6 + 6v 7 + V s . 




(3.9) 



It is now straightforward to obtain the eigenvalues of this two by two matrix. They will not be written here explicitly 
but, as explained above, they are used to calculate the internal energy and specific heat that are displayed in Fig. 

E 

Due to frustration, a finite zero-temperature internal energy is observed for q = 2 in the antiferromagnetic regime 
(v = —1). It is actually not hard to find the exact ground state configuration which consists of only two edges 
per unit cell with energy — J. Since there are five sites per unit cell, the energy per site becomes — 2J/5 = — 0.4J 
as observed in Fig. [8^- For larger values of q there is no frustration and the energy goes to zero as v — > — 1. At 
low temperatures on the ferromagnetic side (v — > oo) the ground state consists always in an alignment of all the 
spins. Thus, each one of the eight edges in the unit cell will have energy — J, resulting in an energy per site of 
—8 J/5 = —1.6 J, as shown in Fig. [5^. 

The specific heat presents a peak in the ferromagnetic side, which is related to the phase transition known to 
be present in the 2D kagome lattice. As the strip widens, the peak will evolve into a characteristic discontinuity of 
the transition. Furthermore, since the Wu conjecture |35j determines with great accuracy that the ferromagnetic 
critical points are located at v c (q = 2) « 1.542, v c (q = 3) w 1.876 and v c (q = 4) « 2.156, we know that the position 
of the peaks will move to the left for wider strips. 

As expected, the Potts antiferromagnet is much more complicated. The q = 2 case on the kagome lattice has 
been solved exactly [39] and it is known to present no phase transition at any temperature. Moreover, it is expected 
that the q = 3 model is critical at zero temperature [40 [ 141] and that it remains noncritical for any q > 3. Thus, 
the peaks for v < observed in Fig. [HJd are not expected to evolve into phase transition divergences as the strip 
widens for q — 2,4; whereas for q = 3 we do expect that to be the case. It is interesting to note that indeed the 
most pronounced maxima corresponds to q = 3. 
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-E/J C/k, 




Figure 8: Results for the kagome lattice strip, (a) Reduced internal energy as a function of the temperature 
parameter v = e K — 1 for different values of q. Lattice frustration gives rise to a finite zero-temperature internal 
energy for q = 2. (b) Dimcnsionlcss specific heat as a function of v. 




Figure 9: A strip of the diced lattice with L y = 2. Each unit cell is composed of five vertices and eight edges. A 
dashed line marks the division between unit cells. 



3.2 Diced lattice strip 

The diced lattice consists on a periodic tiling of the plane by rhombi and is the dual of the kagome lattice. The 
Potts model on the two-dimensional diced lattice has interesting critical properties. For instance, it has been shown 
recently [42] that, even though it admits a height mapping representation for q = 3, it contradicts theoretical 
expectations having q c (diced) > 3. Furthermore, it is the first known case of a two-dimensional bipartite lattice 
with q c > 3. 

We calculate here the exact partition function for the L y = 2 diced lattice strip shown in Fig. [9] Proceeding 
as before, applying the deletion-contraction theorem to all the vertices in a block, the elements of the resulting two 
by two transfer matrix A d (q,v) are, 

af t = q 5 + 8q 4 v + 27q 3 v 2 + 50gV + 52qv 4 + 2q 2 v 4 + 26v 5 + 7qv 5 + 9v 6 + v 7 , (3.10) 

a? 2 = q 4 v 2 +6q 3 v 3 + I6q 2 v 4 + 22qv 5 + q 2 v 5 + I5v 6 + 4qv 6 + 7v 7 + w 8 , (3.11) 

= q 4 + 8q 3 v + 26q 2 v 2 + q a v 2 + Uqv 3 + 6q 2 v 3 + 33w 4 + 19qv 4 + 26v 5 + 2qv 5 + 8v 6 + v 7 , (3.12) 

a\ 2 = q 3 v 2 +6q 2 v 3 + 17qv 4 + q 2 v 4 + 20v 5 + 8qv 5 + I9v 6 + qv 6 + 7v 7 + v 8 . (3.13) 

By the method described above we can now compute the internal energy and specific heat (see Fig. (|10|0. 

In this case there is no frustration and even for q — 2 we observe that the energy reaches E = at zero 
temperature. On the ferromagnetic side, the spins become aligned at zero temperature and the minimal value of 
energy per site E = — 8J/5, is reached. 

For q = 2 the model has been solved exactly in two dimensions and it is known to present an antiferromagnctic 
transition at v ps —0.3401 [39]. Thus, we expect that the observed peak in the specific heat will evolve into the 
corresponding divergence for infinite width. On the other hand, since for q = 3 and 4 the lattice is known to have 
no antiferromagnctic phase transition, we expect that for such values of q the specific heat peaks at v < will 
disappear or move below v — — 1 as the strip widens. 
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Figure 10: Results for the diced lattice strip, (a) Reduced internal energy as a function of the temperature parameter 
v = e K — 1 for different values of q. (b) Dimensionlcss specific heat as a function of v. 

E> 

n 2 1 

Figure 11: A strip of the "shortest path" lattice with L y = 2. Each unit cell is composed of four vertices and eight 
edges. 

3.3 "Shortest-path" lattice strip 

Here we will consider the case of a lattice strip which we will call "shortest path" , since its topology naturally 
appears when one solves the problem of the shortest path connecting the vertices of a square (see Fig. ITTI) . The 
recursive equation method yields the following result for the elements of the transfer matrix A sh (q, v), 

afl = q 4 + 8q 3 v + 27q 2 v 2 + 'i8qv 3 + q 2 v 3 + <i0v 4 + 7qv 4 + 16v 5 + 2v 6 , (3.14) 

afS> = q 3 v 2 + 7q 2 v 3 + 23qv 4 + 35w 5 + 5qv 5 + 26v 6 + 8v 7 + v 8 , (3.15) 

afi = q 3 + ? q 2 v + q 3 v + I8qv 2 + 8q 2 v 2 + 20v 3 + 23qv 3 + q 2 v 3 + 32v 4 + bqv 4 + Uv 5 + 2v 6 , (3.16) 

af 2 = 2q 2 v 2 + I0qv 3 + 2q 2 v 3 + 20v 4 + I3qv 4 + 39v 5 + 3qv 5 + 26i> 6 + 8v 7 + v 8 . (3.17) 

Internal energy and specific heat are plotted in Fig. [12] for different values of q. 

Frustration for q = 2 is manifested in the finite energy at zero temperature in the antiferromagnetic side. As 
for the case of the kagome strip, it is easy to find the minimum energy configuration for this system and find that 
its energy becomes E = — J/4, as can be verified in Fig. [T2Ta). 

Although to our knowledge this lattice has not been studied before, due to universality in the ferromagnetic 
side it is safe to state that the observed specific heat maxima at v > are related to the divergences at the transi- 
tions expected in two dimensions. On the other hand, the Potts antiferromagnet is highly dependent on its lattice 
structure and extracting information about its critical properties would require further study. 

As a further illustration of the power of the present method, in Appendix A we briefly analyze the Fisher zeros 
for the square and the "shortest-path" lattices in the standard range of values of q and also in the large q limit. 




4 Kagome strip L y = 5 

In previous sections we have demonstrated the versatility of the deletion-contraction procedure in handling basically 
any lattice strip of width L y = 2. We will now show that the method can also be applied to wider strips. So far, 
the calculations have been simple enough so that they can be managed by hand. Nevertheless, as the strips become 
wider, the size of the transfer matrices grows rapidly and performing the corresponding calculations requires the 
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assistance of a computer software. In particular, the size of the transfer matrix for a kagome strip of odd width L y 
under free boundary conditions is given by Catalan's number, 

which gives the number of non-crossing partitions of the set {1, 2, £ y }H Thus, for a kagome strip of width L y = 5 
we need to calculate a 42 x 42 transfer matrix. 

To achieve this we have developed a code in Mathematica that is able to handle such a convoluted calculation. 
Basically, the code performs the same task that one would do by hand: Repeatedly apply the deletion contraction 
theorem to one block of the strip. It does so for each one of the 42 possible initial graphs and by keeping track of 
the coefficients corresponding to the final graphs, it constructs the desired transfer matrix. The result is obviously 
too lengthy to display here, but is available upon request from the authors. 

As before, the specific heat can be calculated using (|3.2|) . The result displayed in Fig. Q21 shows some interesting 
features. In the ferromagnetic side the maxima become more pronounced as expected for wider strips. Furthermore, 
all of the peaks come closer to the phase transition points expected for the two dimensional kagome lattice (y c (q = 
2) 1.542, v c (q = 3) « 1.876 and v c (q = 4) « 2.156). 



5 Ornamented L y = 3 hybrid strip 

So far we have dealt with topologically complicated but otherwise regular strips. However, within the same frame- 
work, we can solve even more complex periodic structures such as the configuration depicted in Fig. 1141 It consists 
of a mixture of square and triangular strips decorated by segments of the shortest-path lattice. Obviously, one 
expects that for such a lattice the construction of the exact solution should become much more difficult. There 
are two main reasons for this: Firstly, in this example the lattice has not one but three different types of basic 
boxes (square, triangular and shortest path in the above example). Thus, the system has less symmetries and the 
dimension of the corresponding transfer matrix increases. Secondly, the shape is not regular so that it is not obvious 
how to apply the more usual techniques in this case. Nevertheless, such cases can still be efficiently handled within 
our framework. 

To understand how can one deal with this complicated lattice it is convenient to think of the problem in different 
steps. First, we decompose the ornament of the last block making use of the corresponding transfer matrix calculated 
in this case in Sec. 13.31 It is easy to notice that at the end of this process we end up with the two types of blocks 
shown in Fig. [T47b). Following again the method described in Sec. [2] we can then calculate their corresponding 
5x5 transfer matrices A^> and A^> that enable us to decompose the main body of the strip. Lastly, by combining 
these steps we obtain, 

~t{n) = (a{m)A^ + /3(m)^ 2 >) A^(n - 1), (5.1) 

5 For even widths L y the size is reduced due to reflection symmetry. 
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Figure 13: Specific heat for an infinite kagomc lattice strip of width L y — 5 as a function of v. 
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Figure 14: (a) Right end of a mixture of triangular and square lattices "ornamented" periodically by a series of 
shortest-path strips of length m. (b) Intermediate blocks B^ and Bp) to be considered in the calculation of the 
full transfer matrix of the structure in (a) . 



where a(m) and /3(m) are scalar quantities given by 

a(m) = {(A sh ( (7 , u )) m } ii , p(m) = {(A sh (q,v)) m ) ^ 



(5.2) 



The details of the calculation and explicit expressions for all of the matrices involved can be found in the Appendix. 
Using the same method as in previous cases, we can extract the internal energy and specific heat displayed in Fig. 

□22 

Finally, we point out that the solution outlined here can be straightforwardly generalized for any type of L y = 2 
ornament for which the transfer matrix is known. One only needs to replace a(m) and j3(m) accordingly in (15.11) . 



6 Step-dependent transfer matrix 

In previous sections, it has been shown how to exploit directly the periodic recursive symmetry of lattices. The 
idea so far has been to derive a homogeneous first order linear recurrence relation for a vector which allows to 
solve recursive lattices which are more general than strip lattices. The reason is that our vectorial recursive 
system of equations only need periodic recursive symmetry in order to be applied. Therefore, one can also analyze 
ornamented strip lattices which are, as it has been already explained, strip lattices whose geometry and width 
can change periodically. To achieve this goal, it is very important to analyze all the components of our vectorial 
recursive system of equations, and not just the component which represents the sought partition function. It is 
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Figure 15: Results for the "ornamented" lattice strip in Fig. ITla . (a) Reduced internal energy as a function of the 
temperature parameter v = e K — 1 for different values of q. (b) Dimensionless specific heat as a function of v. 



the interplay between the different components of the vector which provides one with the possibility to build up 
the exact solutions of ornamented strip lattices. In this sense, the present formalism is an extension of the classic 
results of Biggs, et. al. [55] and of Beraha, et. al. [37J [55] m which the authors considered a homogeneous higher 
order linear recurrence relation with "constant coefficients''^ for a scalar in order to find the partition function of 
periodic strip lattices. 

In this section, we will present a system which is quite interesting since it discloses a possibility to further extend 
these results. Up to now, we have only analyzed cases in which the lattice has a periodic recursive symmetry, namely 
the recursive steps to construct the lattice are "always the same" as in the cases of strip lattices as well as ornamented 
strip lattices. Indeed, these are the most analyzed cases in the literature. On the other hand, the key idea of the 
present method also works for recursive non-periodic lattices. We will now show how to apply the present formalism 
to these cases as well. 

An interesting example is the lattice depicted in the Fig. 1161 The resulting linear system will indeed have an 
n-dependent transfer matrix, 

~${n) = A(n)~£(n - 1) (6.1) 

with, 

/ i(T(n + 2) -„»"«) »"« \ 

where T(n) stands for the partition function of the ring with n vertices given by 

n-i 

T(n) = qv" +q^2v l (q + w)"" 1 ^. (6.3) 

i=0 

The n-dependence in (|6.2[) is a manifestation of the fact that each new layer in the lattice is not the same as the 
previous one. Nevertheless, it is possible to apply the method in exactly the same way as before. The vector of the 
partition function can be written in the form 

~t{n) = A(n)A(n - 1)...A(1) if (0), (6.4) 

where the boundary term is given by the two component vector ^(0) = (q,q). The partition function can be read 
from the first component of (|6.4|) . Moreover, it turns to be quite complicated to resolve in general the characteristic 
equation corresponding to the recursive equation of the present example because the coefficients depends on n. 



7 Conclusions and perspectives 

We have exploited a method based on the deletion-contraction theorem to obtain a variety of new exact results for 
the partition function of the Potts model. It has been demonstrated that this technique is quite effective even in 

6 Namely, a linear recurrence relation whose coefficients do not depend explicitly on the recursive step n. 
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Figure 16: Example of a lattice in which the transfer matrix of every layer explicitly depends on the recursive step 
n. 

many cases in which the lattice topology is complicated. Indeed, in order for such a method to be applied, the only 
requirement is that the lattice has a recursive symmetry while there is no need to assume extra symmetries. 

In the case of lattice strips, we have shown that the technique is not restricted to be applied to the simpler 
cases of L y = 2 or 3. With the aid of a computer program we where able to find the exact transfer matrix for the 
L y = 5 kagome strip. In future work we intend to use the same approach to obtain exact results for the L y = 5 
diced lattice strip. 

It has also been shown that this procedure can be easily applied to recursive lattices similar to strips but whose 
width changes periodically (for instance, from Ly = 3 to Ly = m, where m is an arbitrary integer number). The 
possibility to analyze lattices with less regular shapes than strips is a great advantage of the present method and 
could be quite important for applications. 

It is worth noting that this technique can be extended to deal with other boundary conditions such as periodic, 
toroidal or Mobius. Wider strips, which may capture two dimensional features more faithfully, can be considered 
as well. Furthermore, since the calculation procedure is purely combinatoric, it works equally well for non-planar 
graphs. This approach offers also new perspectives for the analysis of scaling and the large q limit. 

Appendix A: Some remarks on Fisher zeros and the large-g limit 

It is sometimes useful to find the zeros of the partition function in the complex v plane. These are the so called 
Fisher zeros which have been used to develop very powerful tools in statistical mechanics [43l |44j H3 [46 , 47, 48 . 
We will show the Fisher zeros of the square and "shortest-path" strip lattices for two ranges of values of q: small 
integers and the large q limit. 

When q is in the small integer range the structure of the Fisher zeros for the two lattices is very sensible to the 
topology of the lattice as illustrated in Fig. [17] On the other hand, when q is very large it is worth noting that the 
curve in the complex plane simplifies to a circle (see Fig. I18[) . 

Let us now consider the recursive lattice studied in Sec. [6] The n-dependent transfer matrix depends explicitly 
on T(n) which is the partition function of a ring made of n segments. Therefore, it is natural to expect that 
the specific heat of this lattice will present many peaks corresponding to the various factors T(n') (n' = 1, ...,n) 
appearing in the exact expression of the partition function. This expectation is confirmed in Fig. 1191 where the 
large-g specific heat displays these features. 

Appendix B: Ornamented strip transfer matrix 

As explained in the main text, one can proceed by first decomposing the ornament and then the main body of the 
strip. Combining both of these procedures the resulting transfer matrix is given by, 



where a(m) and /3(m) are scalar quantities that can be extracted directly from the transfer matrix of the "orna- 
mented" strip, 




(7.1) 




(7.2) 
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Figure 17: The Fisher zeros are plotted in the complex v plane for strips of length n = 100 for the square lattice in 
the cases (a) q = 2 (b) q = 3 and (c) q = 5; and for the "shortest-path" in cases (d) q = 2, (e) q = 3 and (f) g = 5. 
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Figure 18: The Fisher zeros are plotted for strips of length n = 100 and q = 2 200 for: (a) square lattice and (b) 
"shortest-path" lattice. 



On the other hand, A^> and are 5x5 matrices given by 



OCl,C2,C3 
S1,S 2 


r<ci,C3 
^s 2 


nc 2 ,C4,C5 
S1,S 2 


pc 2 ,C4,c 5 
s 2 


FC2 

s 2 


r>c 1 ,c 2 ,c 3 

B t!,t 2 


/^Ci,C3 


7-)C2,C4,C5 


A t2 


^*2 


OcZi ,d 2 ,<^3 
«1,«2 


£>di,d 3 


r\d 2 ,d4,d s 

Sl,S 2 


S<2 


pd 2 

S 2 


B t u t 2 


/~idi,d 3 

°*2 


j~ ) d 2 ,d i ,d !i 
U t^,t 2 


A *2 


pd 2 
r t 2 


,e 2 ,e 3 ,e 4 

,c 2 ,C3',di ,d 2 : d$ 


Ce 3 ,e 4 


fjei,e 2 

c 2 ,C4,c 5 ;d 2 ,d4,d 5 


J -'e3,e4 


P 

r e 3 ,e4 



A™ = (l + v) 



/ nci,c 2 ,c 3 

' S3,S4 


s 4 


nc 2 ,C4,C5 
S 3 ,S4 


PC2,C4,C 5 
S4 


PC2 

S4 




-gci,c 2 ,c 3 

^3,^4 




^)C2,C4,C5 
*3,*4 


77iC 2 ,C4,C5 
*4 


PC2 
r «4 




r>d 1: d 2 ,d 3 


/"idi,d 3 


f)d 2 ,d4,d 5 


p°d 2 ,d4,d 5 


pd 2 




S 3 ,S4 


S4 


S3.S4 


S4 


S4 




r>d 1: d 2 ,d 3 

D H,t4 


°*4 


r^d 2 ,d4,d$ 

U t 3 ,t4 


A *4 


pd 2 




i -g/1,/2,/3,/4 

\ Ci,C 2 ,C3,di,d2,c!3 


C h,f4 


fS/li/2 

C2,C4,C5,d2,d 4 ,(i5 


E hj4 


F h,h 


) 



(7.3) 



(7.4) 
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Figure 19: (a) Specific heat in the ferromagnetic range for the n = 50 lattice. A global maximum and many local 
peaks in the high temperature regime are observed, (b) Lower temperature regime in more detail. The density of 
local peaks is observed to increase with temperature. 



where. 



DCl,C2,C3 

si,s 2 

r)C2,C4,C5 
Sl,S2 

s 2 

S2 

Sei,e2,e3,e4 

Cl,C2,C3,d 1 ,d 2 ,d,3 



2^ei,e2 

C2,C4,C5, £02,(^4, ^5 



with T 2 = 9(9 + 2i>) + gw 2 , and 



si((q + u)ci + c 2 + (1 + «)c 3 ) + s 2 Ci, 

s 2 (vci + (1 + u)c 3 ), 

si(vc 2 + (q + v)c 4 + (1 + v)c 5 ) + s 2 c 4 , 

s 2 (vc 2 + vc 4 + (1 + v)c 5 ), 

S2C2, 

ei((g + u)ci + c 2 + (1 + v)c 3 ) + e 2 ((q + v)d t + rf 2 + (1 + u)of 3 ) 

+(g + 2v)e 3 + (1 + v)ei, 

v(v + T 2 /q)e 3 + v(l + v)(2 + v)e 4 , 

ei(fc 2 + (q + v)c 4 + (1 + v)c 5 ) + e 2 {vd 2 + (q + v)d 4 + (1 + v)d 5 ), 
v 2 (v + T 2 /q)e 3 + v 2 (l + v)(2 + v)e 4 , 
v(q + 2v)e 3 + v(l + v)e 4 , 



si = q + 2v, s 2 =v 2 , s 3 = 1, 



,S 4 = V, 



h = l + v, t2=v(l+v), t 3 — 0, t 4 = l + v, 
c 4 =q + 3v, c 2 =v 2 , C3=v 2 , c 4 = v 2 , C5 = v 3 , 
di = (l + v), d 2 =v(l + v), d 3 =v(l + v), d 4 = 0, d 5 = v 2 {l + v), 
e\ = 1, e 2 = v, e 3 = v, e 4 = v 2 , 
fi = 0, /a = 0, / 3 = 1, /* = »• 

As usual the general expression for the vector of partition functions is given by 



~t(n) = [(a(m),4 (1) +£(m)yl (2) ) A^j" ^(1) 



(7.5) 
(7.6) 
(7.7) 
(7.8) 
(7.9) 
(7.10) 

(7.11) 
(7.12) 
(7.13) 
(7.14) 
(7.15) 



(7.16) 
(7.17) 
(7.18) 
(7.19) 
(7.20) 
(7.21) 
(7.22) 



(7.23) 
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where the initial data correspond to the vector 



/ 



q(q + v) 2 \ 
q{q + v)(l + v) 
q{q + v){l + v) 



1(1) 



(7.24) 



V 1 



q{l + v) 2 
(q + v) + qv(l + v) ) 
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